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Abstract 

In the standard framework of non-equilibrium quantum field theories, 
the pinch singularities associated to multiple products of ^-functions do not 
cancel in a perturbative expansion unless the particle distributions are those 
for a system in thermal and chemical equilibrium. 
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Since the pioneering work of Keldysh |T|, there have been many attempts 
to develop a formalism for non-equilibrium systems. Most of the studies 
follow the original method of Schwinger 0, which uses a Closed-Time-Path 
(CTP) in the complex-time plane |§. Unlike the imaginary-time formalism, 
but closely related to the real-time formalism, this approach leads to a 2 x 2 
matrix stucture for the propagator 

One other approach, which is more recent, is Thermo-Field Dynamics |J. 
There, the doubling of the degrees of freedom is assumed from the beginning. 
For equilibrium systems, one can show that the theory is completely equiva- 
lent to the real-time formalism Q . The non-equilibrium version of TFD has 
been intensively studied over the last few years ||. 

In this letter, we shall not discuss the relative advantages of one method 
over the other, but shall just make a very simple remark about the pertur- 
bative expansion, when one is using non-equilibrium particle distributions. 
We consider the scalar case. The most common form of the free propagator 
that is used in the literature is |3|, |J 

f D 11 (K) D 12 {K) 
v D 21 (K) D 22 (K) 




+ 



n(Jfe) 6(k )+n(k) 
y9(-k )+n(k) n(k) 

with the usual vacuum Feynman propagator 



2n8(K 2 - mil) 



A ( K ) = K 2 (2) 
K z — m z + ze 

and where n(k) is an arbitrary distribution function (positive-definite). As a 
matter of fact, this is the only difference with the equilibrium case, where the 
propagator is exactly the same as in eq. (|1|), except that n(k) = l/(e' fco '/ T — 1). 



Self-interactions are taken as in the equilibrium case, i.e. one needs to 
distinguish two types of vertices, which will be called type-1 and type-2 ver- 
tices in the following. For a Feynman diagram with a certain configuration of 
type-1 and type-2 vertices, one uses a D a b(K) propagator when the momen- 
tum K flows from a type-a vertex to a type-6 one. For a g n cj) n fn\ interaction, 
type-1 vertices get a factor {ig), whereas type-2 vertices get an opposite fac- 
tor, (— ig). These Feynman rules set up a formalism for non-equilibrium 
systems. They have been used many times in the literature [1-9]. 

Let us now see if the above Feynman rules give a well-defined perturba- 
tive expansion. Knowing the problems which arise at equilibrium, the first 
obvious question is the absence of pathologies, or "pinch singularities" . This 
is due to the presence of several 5(K 2 — m 2 ) terms in eq. (1). They appear 
as multiple products in high order calculations, and for individual graphs, 
lead to mathematically ill-defined expressions. At equilibrium, one has to 
sum over the different types of vertices in order to obtain the cancellation of 
such pathologies [|j]. Surprisingly, this very question has almost always been 
occulted for the non-equilibrium case JF], ||. 

We first consider the simplest case where such a cancellation should occur, 
if any. Suppose that we want to calculate the tadpole contribution at the 
2-loop level shown in fig. 1, in the g(</> 4 )4 interaction model. It is given by 

° 2 1 I l^f [{Dn{P)?D u {K) - D 12 (P)D 21 (P)D 22 (K)} . (3) 

Observing that the tadpole does not have any imaginary part and that 
KeDu(K) = KeD 22 (K), one can factorize the i^-integral and obtain 

g2 I !0f KeDn{ - K) I [(^( P )) 2 " D 12 (P)D 21 (P)} . (4) 



Then an easy algebraic manipulation shows that the term in the square 
brackets simplifies to 

[...] = (A(P)) 2 + n(p) ((A(P)) 2 - (A*(P)) 2 ) , (5) 

which shows the absence of pinch singularities, or A(P)A*(P) terms in the 
final result. The cancellation procedure works essentially the same way as 
in the equilibrium case and is here totally independent of the distribution 
function n(p). This result had previously been found in J7J. 

One should not claim victory too soon, though. This simple exercise in 
the (p 4 model does not illustrate well the game of cancellation that is at play. 
In particular, in the previous example only two terms are involved, although 
the general case involves four terms. 

Let us next consider a more complicated case. The only place where the 
pinch singularities can occur is in repeated self-energy insertions. Consider 
for instance the diagram shown in fig. 2. It contains the following expression 

Y,D la (P)Z ab (P)D b2 (P) = P 11 (P)S 11 (P)P 12 (P) +P n (P)S 12 (P)P 22 (P) 

+P 12 (P)S 21 (P)P 12 (P) +P 12 (P)S 22 (P)P 22 (P), (6) 

which must be free of pinch singularities, as it enters directly, at the two-loop 
level, into the calculation of a physical quantity (the decay rate), as we shall 
see in the following. For the same reason, the cancellation must also take 
place separately for p > and for p < 0. 

The different components of the self-energy can be related to each other 

by 

S n (P) = -E*,(P) 
ImSn(P) = ^(S 12 (P) + S 21 (P)). (7) 
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These relations follow from the definition of the two-point matrix Green's 
function in different chronological products, using the standard CTP contour 
|||. They are independent of perturbation theory. Then, using (0), one can 
show that eq. (H) is free of pinch singularities provided 

(6{p )n(p) - e(-p ){l + n{p))) £ 12 (P) = e(p )(8{p ) + n(p))£ 21 (P). (8) 

For p > 0, one has 

n(p)£ 12 (P) = (l + n(p))£ 21 (P), (9) 

and when n(p) = l/(e P0//T — 1), this gives 

£ 12 (P) = e^ T £ 21 (P), (10) 

which shows that eq. (H) can in fact be regarded as a non-equilibrium exten- 
sion of the KMS relation. 

The quantities S i2 and £ 2 i, which are related through eq. (||), are two 
independent physical quantities. For a A0 3 interaction model in n space-time 
dimensions, one has, at one loop, 

-iS 12 (P) = X 2 Jj^- 2 (e(k )+n(k))(e(p -k )+n(p-k)) 

5(K 2 — m 2 )5((P — K) 2 — m 2 ). (11) 

The kinematics are the same as in the equilibrium case, i.e. (E —p)/2 < 
ko < {E + p)/2, for po = E > and P 2 > m 2 . Then, the statistical factors 
in the above equation are just the ones corresponding to outgoing particles. 
As in the equilibrium case, this allows us to relate Si 2 with the absorption 
(or decay) rate of the particle 0. Similarly 

r d n K 

-zE 21 (P) = A 2 J ^^-^(e(-k ) + n(k))(9(-po + k ) + n(p - k)) 

5(K 2 — m 2 )5((P — K) 2 — m 2 ) (12) 



is related to the emission (or creation) rate. Strictly speaking, when P is 
on shell, there is no kinematical phase space anymore in eqs. (|TT|) and (|T^) 
and both expressions are equal to zero. But this is only true at the one-loop 
level. Beyond one-loop calculations, £12 and £21 no longer vanish when the 
external momentum goes on shell. Another way of looking at the problem is 
to consider two types of particles, <p\ and 02, with two different masses and 
a coupling such as X(f)i4>2 (f° r details see ||). In this case, the interpreta- 
tion of £12 and £21 as, respectively, the decay and emission rates, is more 
transparent. 

According to the above discussion, the time evolution of the particle num- 
ber density follows 

- 2zpo^|^ = (1 + n(p))J^(P) - n{p)H 12 {P). (13) 

We see that in order to have the cancellation of pinch singularities, one must 
have 

= 0, (14) 

at 

which is quite disappointing for a non-equilibrium framework ! This is the 
first contradiction. However, this is not a real problem as it can be realized 
that for the propagator defined in (1), it must be assumed that the time 
variation of the density matrix is slow compared with the typical time scale 
between the particle interactions. If this is not the case, the Fourier trans- 
form in (1) just does not make sense. In the context of using (1), there is 
reversibility at the microscopic level, but one can impose some irreversibility 
at the macroscopic level. 

But the main trouble is that, even assuming a slow variation of the density 
matrix, eq. ([|) is not guaranteed to hold. The point is that the micro- 
reversibility conditions are well-known to be satisfied only by equilibrium 
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distributions [K|. This can be verified by explicit calculations using eqs. 
(JTTJ) and (|l^). Equation ^) can only be satisfied if the distribution n(p) is 
of the Bose-Einstein type. The only alternative to this result is to give up 
energy conservation at the vertices, which is not very satisfactory. 

To see how deeply rooted the problem is, let us consider again the case 
of two types of particles, with two different initial temperatures T\ and T2, 
and a single weak interaction A0102, which is switched on at some arbitrary 
time to- Then the free propagators for 0! and 02 are just the same as in 
eq. (1), with n(p) the Bose-Einstein distribution, but different temperatures 
for 0i and 02- At one loop, the self-energy for 0i involves only 02 fields, and 
it obeys the relation 

S 12 (P) = e^S 21 (P), (15) 

as 02 is thermalized with temperature T 2 . On the other hand, the temper- 
ature that enters n{p) in eq. (H) is 7\, not T 2 , so that the cancellation of 
pathologies occurs only when 7\ = T 2 . 

Perhaps the most disappointing lesson of all this is to realize that it is not 
even possible to look at small deviations from equilibrium. Also, this problem 
is not specific to the relativistic case: it shows up in the same fashion in the 
non-relativistic limit. 

Except for |§, this fact seems to have gone unnoticed in the literature. 
One way of solving this problem is to use Schwinger-Dyson equations. By 
using only exact propagators in one-loop calculations, there is no possibility 
of having pinch singularities. This is equivalent to giving up perturbation 
theory. As a matter of fact, this has been the usual way of doing calculations 
for non-equilibrium systems [[|. However, at some point, and unless the 
problem can be solved exactly, which is rather rare, one is forced to use some 
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perturbative input. In the light of our results, this has very weak justification 
as the bare perturbation series is ill-defined. In particular, terms that are 
associated with pinch singularities (which can be regularized by introducing 
some finite width in the propagator) are likely to give large and uncontrollable 
contributions. 

In conclusion, it is impossible to make use of perturbation series with 
the propagators (1) outside an equilibrium framework. One must use time- 
independent Bose-Einstein distribution functions (the same is true for fermion 
fields, which have to obey the Fermi-Dirac statistics). This guarantees the 
cancellation of pinch singularities at all orders of the perturbation series. This 
cancellation is intimately tied to the micro-reversibility conditions. Note also 
that the same conditions ensure the cancellation of infrared and mass sin- 
gularities (KLN theorem), which defines a well-behaved perturbation series 

0- 

The problem can clearly be solved in principle by working in the T = 
representation of the system, but then the calculation is computationally in- 
tractable because of the complexity of the system. If we choose to work with 
any fixed-T representation, then the system quickly leaves the vacuum state. 
We cannot treat the problem using the propagators (1), because then the 
5-function pathologies discussed here arise, so our only choice is to allow the 
state of the system to depart from the thermal vacuum; again, the problem 
quickly becomes computationally intractable because of the complexity of 
the state. Finally, we could choose the closest finite-T vacuum to approxi- 
mate the state of the system at each time, but then we are faced with the 
intractable task of transforming the state of the system as a function of time 
as the Fock states change with the changing thermal vacua. 
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We stress again that anything perturbative beyond linear response theory 
does not seem easily realizable. At present, there does not exist any correct 
way for deriving a consistent perturbative Green's function formalism for 
a system that is even slightly out of equilibrium, without losing all of the 
advantages of working in finite-T vacua. 
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Figure captions 

Fig. 1 A 2-loop contribution to the self-energy. 

Fig. 2 A particular summation of self-energy insertions. 
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This figure "figl-l.png" is available in "png" format from: 
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